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1.  In  a paper  published  under  the  above  title  in  the  Philosophical 
Magazine  and  Journal  of  Science  for  1892,  Lord  Kelvin  described  a 
graphic  method  for  the  solution  of  dynamical  problems.  We  will  quote 
from  his  paper  the  following  account  of  his  method. 

«In  dynamical  problems  regarding  the  motion  of  a single  particle  in 
a plane  the  following  plan  is  given  for  drawing  any  possible  path  under 
the  influence  of  a force  of  which  the  potential  is  given  for  every  point  of 
the  plane.  Suppose,  for  example,  it  is  required  to  find  the  path  of  a par- 
ticle projected,  with  any  given  velocity,  in  any  given  direction  through 
any  given  point  P{)  (see  fig.  i)».  L 

« Calculate  the  normal  component  force  at  this 
point,  and  divide  the  square  of  the  velocity  by  this 
value,  to  find  the  radius  of  curvature  of  the  path  at 
that  point.  Taking  this  radius  on  the  compasses, 
find  the  centre  of  curvature  C0  in  the  line  PqK,  per- 
pendicular to  the  given  direction  through  P0,  and  de- 
scribe a small  arc,  P^PXQX,  making  PxQi  equal  to 
about  half  the  length  intended  for  the  second  arc. 

Calculate  the  altered  velocity  for  the  position  0 x, 
according  to  the  potential  law,  and,  as  before  for  Pot 
calculate  a fresh  radius  of  curvature  for  Ql  by  find- 
ing the  normal  component  force  for  the  altered  direc- 
tion of  normal,  and  for  the  velocity  corresponding  to 
the  position  of  Q L.  With  this  radius,  find  the  po- 
sition of  the  centre  of  curvature,  Ct,  in  P^C^L,  the 
line  of  the  radius  through  P\.  With  this  centre  of 
curvature  and  the  fresh  radius  of  curvature,  describe  an  arc  PXP2Q2) 
making  P2Q 2 equal  to  about  half  the  length  intended  for  the  third  arc; 
calculate  radius  of  curvature  for  position  Q.,,  draw  an  arc  P^P^Qz*  and 
continue*  the  procedure ». 
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The  practical  utility  of  this  method  is  dependent  upon  the  ease  or 
difficulty  with  which  the  radius  of  curvature  can  be  calculated.  If  the 
potential  on  which  the  force  depends  is  a somewhat  complicated  function 
of  the  coordinates  of  the  point,  the  method  immediately  becomes  so 
troublesome  to  use,  that  recourse  should  rather  be  had  to  one  of  the  many 
numerical  methods  of  integration  for  the  calculation  of  the  path. 

We  intend,  however,  to  demonstrate  an  extremely  simple,  approximate 
construction  of  the  radius  of  curvature  by  the  aid  of  a number  of  cquipo- 
tential  lines,  whereby  the  above  difficulty  is  avoided. 

Let  the  differential  equations  for  the  motion  of  the  point  be 


d2x  _dU 
dt2  ~ dx 


d2y_dJJ 
dt2  dy 


(I) 


where  U is  a given  function  of  x and  y.  Let  us  further  assume  that 
we  have  drawn  in  advance  a number  of  equipotential  lines, 


U=0, 


U = 


Let  A M B be  a part  of  the  path,  and  M its  point  of  intersection 

with  the  equipotential  line  £/=— . Further  let  M N be  the  normal,  with 

n 

direction  towards  increasing  O,  i.  e.  towards  the  centre  of  curvature  of 
the  path.  If  we  designate  as  D the  part  of  the  normal  between  the  two 

nearest  equipotential  lines-,  U = - — — 1 and  U — we  b(lve  for  the 

radius  of  curvature  q in  the  point  M,  the  approximate  formula 


PD . 


By  the  aid  of  this  formula,  the  radius  of  curvature  may  be  easily  and 
quickly  constructed,  when  the  equipotential  lines  have  been  drawn  once  for 
all.  The  construction  of  the  path  of  the  point  by  the  aid  of  the  radii  of 
curvature  thereby  becomes  very  simple  and  easy  of  accomplishment.  We 
will  return  later  to  the  question  of  the  possibility  of  modifying  the  method 
when  the  radius  of  curvature  is  very  large,  or  when  the  equipotential 
lines  lie  far  from  one  another. 
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We  will  first  prove  the  correctness  of  the  approximate  formula  q — pD. 
If  we  designate  the  projection  on  the  normal  M N of  the  force  acting  on 
the  point  M,  as  Kn}  and  the  velocity  of  the  point  as  v,  we  have,  as  we 
know,  the  mass  of  the  point  being  equal  to  1, 


dU  m 

According  to  system  I,  v2  = 2 Uf  and  Kn  — — y,  where  — , indicates  the 


d_U 

IN 


pf 


derivative  of  U along  the  normal  M N.  This  gives 

1 1 dU 

q 2 U 2N 


(II) 


Let  us  now  consider  an  assumed  point  Mf  on  the  normal,  and  call 
its  distance  from  a fixed  point  on  the  latter  5,  J to  be  reckoned  positive 
in  the  direction  M N.  If  M ' moves,  U will  be  a function  of  5,  and  if 
this  function  be  developed  in  series,  we  easily  obtain  the  desired  approxi- 
mate formula,  o = pD.  It  was  in  this  way,  moreover,  that  I first  demon- 
strated this  formula. 
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We  attain  our  end  most  rapidly,  however,  in  accordance  with  Pro- 
fessor Runge’s  suggestion,  by  regarding  J as  a function  of  U,  and  con- 
sidering the  corresponding  curve, 

i=f(U) 


where  U is  an  abscissa  and  £ an  ordinate  l. 

three  points  whose  abscissae  are  ^ ^ and 

n n 

parallel  with  the  axis  of  abscissae;  then  A A'  = 


Let  A}  B and  C be  the 

^ , and  let  A A'  be 

n 

2,  and  A'  C = D. 
n 


Thus  if,  instead  of  the  direction  of  the  tangent  in  the  point  B,  we 
take  as  approximate  direction  A C,  we  obtain 

. n „ 

tan  ft  — - D, 


ft  being  the  angle  made  by  the  tangent  with  the  axis  of  abscissae. 

On  the  other  hand,  tan  ft  — f (U),  where  the  value  of  U is  — . 
Now  since 


r (U)  = 


dU 

dN 


I am  indebted  to  Professor  Runge  for  information  respecting  this  method  of  procedure. 
At  the  beginning  of  March,  1907,  I had  sent  a description  of  the  method  with  proofs 
to  Professor  Mittag  Leffler,  asking  if  it  were  known.  He  sent  it  to  Professor 
Runge  in  Gottingen,  who  drew  attention  to  the  above  simpler  proof. 
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we  also  obtain  for  tan  ft  the  value  If  we  compare  these  two  expres- 

sions for  tan  [S,  we  obtain 


q = pD 


Q.  E.D. 


It  is  easy  to  find  an  expression  for  the  accuracy  attained,  for,  accord- 
ing to  Taylor’s  formula,  we  have 

f(U,  + JU)  = f(UJ  + f'(U0)jU+fC~^  UU)*  + /'"  - V U)3 


i x 2 


1 X 2 X 3 


and 


/•(£/„  - JU)  = f(Ut)  - f (JUf  - r (JU)S, 

where  6 and  O'  are  numbers  between  o and  1.  From  this  we  get 

f(ut+ju)-f(ut-ju)  — *f’(utuu+wm{ut+e.ju)+r{ut-tf.<iu)\  (jU)3. 

If  we  now  choose  U0  — ^ and  J U — , we  obtain 


D= 0 + a 

p ' 6nd 


+r\f-;) 


and  thus 


where  we  have  put 


Q = pD  + e, 


Now  if  f"(U)  in  the  interval  — — - U I,  lies  between  — m 

n n 

and  -|-  m,  then 

— -^5  m <C  « <C  m • 

3 rr  3 n6 

t a - n a ( t t\  u , , dU  d2U  _nH  d*U 

In  order  to  find  m,  f [U)  may  be  expressed  by  ana  ^3  • 

As  | = f (U),  we  have 


f"  (U)  = 


Yrf8v\2 

W$2/  rf|  rfi 


rfs£/ 

3 


^y\5 
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Thus  the  first  three  derivatives  of  the  function  U along  the  normal  M N 
may  be  formed,  and  (U)  may  then  be  found  by  substitution  in  the 
above  formula. 

We  see  that  the  error  e diminishes  very  rapidly  with  increasing  n. 

2.  The  above  method  was  worked  out  for  the  purpose  of  simpli- 
fying the  study  of  the  integral  curves  of  the  system 

d2R  _ 1 dQ 
ds1  ~*dR 


where  Q is  the  function 

n [2  y R l2 

^ 1 |tf  +«*)*. 

and  y is  constant. 

One  is  led  to  this  system,  when  it  is  a question  of  calculating  the 
trajectories  of  electric  corpuscles  in  space  under  the  influence  of  terrestrial 
magnetism  E 

One  of  my  assistants,  Mr.  E.  Steenstrup,  engineer,  has  constructed 
under  my  guidance,  after  the  above  method,  more  than  50  paths  of  system 
(III),  answering  to  various  values  of  the  constant  y.  The  practical  utility 
of  the  method  has  thereby  been  sufficiently  tested.  The  results  of  these 
constructions  will  be  fully  described  in  Professor  Birk eland’s  report  of 
the  Norwegian  auroral  expedition  of  1902 — 3,  in  the  section  that  treats  of 
the  mathematical  theory  of  auroras  and  magnetic  perturbations. 

3.  We  will  here  briefly  relate  our  experience  with  regard  to  the 
employment  of  our  graphic  method  of  finding  the  paths  defined  by 
system  (I). 

First  a number  of  equipotential  curves, 

U = o,  U = \ 11=-,  U = 3-, 

n n n 


See  my  paper:  „Sur  les  trajectoires  des  eorpuscules  electrises  dans  l’espace  sous  Taction 
du  magnetisme  terrestre,  avec  application  aux  aurores  boreales".  Archives  des  Sciences 
Physiques  et  Naturelles;  July — October,  1907.  Geneva. 
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are  calculated,  and  the  curves  drawn  upon  paper,  n is  assumed  so  large 
that  the  curves  on  the  diagram  are  not  too  far  distant  from  one  another; 
in  some  places  it  will  generally  be  necessary  to  assume  larger  values  for 
n than  at  other  places.  For  the  sake  of  clearness,  it  is  recommended  that 

all  curves  U = - be  dotted  lines  for  p odd,  and  continuous  lines  for  p 
n 

even,  and  that  the  odd  numbers,  1,  3,  5,  7,  etc.,  be  written  in  the  spaces 
between  the  continuous  curves,  from  U = o onwards,  as  shown  in  figure  4. 


If  the  point  M of  the  path  falls  upon  one  of  the  dotted  curves,  the 
radius  of  curvature  then  equals  the  portion  D of  the  normal  between  the 
two  nearest  continuous  curves  multiplied  by  the  number  of  the  interval. 
If  M falls  upon  one  of  the  continuous  curves,  we  have  to  multiply  D by 
the  even  number  lying  between  the  numbers  of  the  adjacent  intervals. 

The  graphic  construction  of  the  path  may  now  be  accomplished  in 
the  manner  given  by  Lord  Kelvin;  but  in  addition  to  compasses,  it  is  an 
advantage  to  use  a narrow  tape  divided  into  millimetres,  by  which  the 
centre  of  curvature  can  be  quickly  marked.  By  means  of  pins,  this  tape 
is  attached  to,  and  detached  from,  the  points  in  the  path  and  their  centres 
of  curvature  alternately,  as  the  construction  progresses.  The  product 
pD  may  be  found  by  the  aid  of  a sliding-rule. 

More  accurate  results  are  obtained  by  substituting  transparent  gelatine- 
paper,  cross-ruled,  for  the  tape. 

Before  proceeding  to  a treatment  of  the  cases  in  which  the  method 
must  be  modified,  we  will  consider  the  connection  between  the  radius  of 
curvature  q of  an  assumed  path  through  the  point  M (see  fig.  5),  and  the 
radius  of  curvature  £0  of  a path  through  the  same  point,  tangent  to  the 
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equipotential  line  through  that  point.  If  we  call  the  projection  of  the 
force  K along  M N,  as  before,  Kn,  we  have 


Now  since 


and 


we  obtain 


Kn  — K cos  ip  . 

K = = 2U 

Q o Q o 


Q 


= -lo- 
cos Ip  } 


whence  we  have  the  construction  of  q shown  in  fig.  5,  when  o0  and  the 
angle  xp  are  given. 


This  construction  is-  to  be  recommended,  when  there  is  a suspicion 
that  the  formula  q = pD  does  not  give  exact  results,  owing  to  the  fact  that 
the  path  is  becoming  perpendicular  to  the  equipotential  curves , or  because 
the  course  of  these  curves  is  not  sufficiently  regular  in  the  vicinity  of  the 
point  M.  By  the  aid  of  the  already-mentioned  cross-rule  gelatine-paper, 

the  factor  - is  quickly  found  as  follows, 
cos  ip 

A cross  of  axes  is  made  on  the  gelatine-paper,  consisting  of  two 
axes  at  right  angles  to  one  another,  one  an  axis  of  abscissae  xx"f  and 
the  other  an  axis  of  ordinates  y y' . A needle  is  then  put  through  the 
intersection  of  the  axes  and  the  point  M in  the  path,  and  the  gelatine- 
paper  is  turned  in  such  a manner  that  the  axis  xx"  is  tangent  to  the 
path.  Along  the  normal  M N,  a portion,  M N’,  is  then  marked  off,  of  a 
given  length  a,  and  N*  is  marked  upon  the  paper  by  making  a hole  with 
a needle  through  the  gelatine.  The  gelatine  is  then  turned  about  M until 
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the  axis  -xx"  is  tangent  to  the  equipotential  curve  through  M,  and  the 
distance,  b,  from  the  point  marked  on  the  paper  to  the  axis  x x' , is  read 


off.  Then  cos  lb  — and  thus 
r a 


Q 0 


whereupon  q is  quickly  calculated  with  a sliding-rule. 

In  the  case  of  q being  so  large  that  the  centre  of  curvature  falls  out- 
side the  paper,  the  following  method  may  be  employed. 


t * 


‘-T-'.-w 


Fig.  6. 

In  this  figure,  where  M M'  is  an  arc  of  a circle  with  radius  q . and 
centre  in  O (centre  of  curvature),  we  have 

rj  = q (1  — cos  a) 
c = q sin  « 

consequently,  with  sufficient  accuracy, 


2 Q 


1 a 
h — — c 


Here  c is  assumed,  as  also  a.  The  position  of  M'  is  thereby  determined, 
rj  and  h being  easily  calculated  with  a sliding-rule.  As  also  the  position 
of  P,  is  known,  we  have  the  direction  of  the  radius  of  curvature  in  the 
point  Mr. 

If  the  gelatine-paper  is  used,  it  is  turned  about  M so  that  the  one 
axis  falls  along  M N' . By  the  aid  of  the  cross-rule  paper,  the  points  M' 
and  Pf  may  then  easily  be  marked. 

This  is  the  description,  in  its  principal  details,  of  the  graphic  con- 
strqction  of  the  path. 
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